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Abstract. We report on observations of SN 1987A with 
ISO SWS/LWS made 9-11 years after the explosion. 
No emission from the supernova was seen. In particu- 
lar, the upper limits on the fluxes of [Fe I] 24.05/im and 
[Fe II] 25.99/im on day 3 999 are - l.lJy and 1.4 Jy, re- 
spectively. Assuming a homogeneous distribution of '^^Ti 
inside 2 000 km s~^, we have made theoretical models 
to estimate the mass of ejected '^'^Ti. Assessing various 
uncertainties of the model, we obtain an upper limit of 
~ 1.5 X 10^^ M0. The implications of this are discussed. 

The LWS data display continuum emission as well as 
nebular lines of [O I], [C II] and [O III] from neighboring 
photoexcited regions in the LMC. The [O III] lines indicate 
an electron density of 120 ± 75 cm~^, and the continuum 
can be explained by dust with a temperature of ^ 37 K. A 
second dust component with ^ 10 K may also be present. 



Key words: supernova: individual: SN 1987A - nucle- 
osynthesis - supernovae: general 



1. Introduction 

Supernova (SN) 1987A provides a convenient tool to test 
our understanding of nucleosynthesis in massive stars and 
during supernova explosions (Thielemann et al. 1996). In 
particular, the late {t ^ 150 days) emission probes directly 
the elemental abundances deep in the stellar ejecta. As 
in other core-collapse SNe without a dense circumstellar 
medium, the energy production is at this epoch dominated 
by radioactive energy from the decay of ^^Co to ^^Fe, the 
cobalt itself being the decay product of ^^Ni. It was quickly 
realized that the ejected mass of ^^Ni in SN 1987A was at 
least 0.05 M0 (Woosley et al. 1988), the actual number 

* ISO is an ESA project with instruments funded by ESA 
Member States (especially the PI countries: France, Germany, 
the Netherlands and the United Kingdom) and with the par- 
ticipation of ISAS and NASA. 
Correspondence to: peter@astro.su.se 



being close to 0.07 Mq (e.g., Suntzeff & Bouchet 1990). 
This decay continued to power the bolometric light curve 
in an undisputed fashion for about ~ 800 — 1 000 days 
(Bouchet et al. 1996). 

However, during the following epoch, radioactive de- 
cay seemed to be unable to explain the bolometric flux 
of the supernova, unless a large amount of ^^Ni (decaying 
to '^^Co and further to ^^Fe) was included (Suntzeff et al. 
1992). Fransson & Kozma (1993) demonstrated that the 
derived ^^Ni/'^^Ni ratio could be much closer to the solar 
^^Fe/^^Fe ratio when time dependence was accounted for. 
This effect, now known as the "freeze-out" effect, stems 
from the fact that the energy stored at earlier epochs in 
the low-density hydrogen gas slowly leaks out and domi- 
nates the bolometric light curve after ^ 1 000 days (Frans- 
son & Kozma 1999). The freeze-out effect becomes less 
important as the radioactive isotope ^''Ti starts to domi- 
nate. This occurs after ~ 1 500 — 2 000 days (Woosley et 
al. 1989; Kumagai et al. 1991; Fransson & Kozma 1999). 
''^Ti decays to ^''Sc on a time scale of 87.0 ±1.9 years 
(Ahmad et al. 1998; Gorres et al. 1998), and then quickly 
further to ^"'Ca. 

It is important to determine the ^^Ni, '^'^Ni and '^'^Ti 
masses in order to constrain models for the supernova ex- 
plosion and the explosive nucleosynthesis (e.g., Timmes 
et al. 1996). In the case of ^''Ti, models predict that 
~ 10"'* Mq could be synthesized in SN 1987A (e.g., Ku- 
magai et al. 1991; Woosley & Hoffman 1991). 

Chugai et al. (1997) estimate that the mass of ejected 
^''Ti should be M(44Ti) ~ (1 - 2) x lO"'* M© from the op- 
tical line emission at 2875 days. Kozma (1999) and Kozma 
& Fransson (private communication 1999; hereafter KF99) 
preliminary find that MC^^Ti) = (1.5 ± 1.0) x lO"'' M© 
best explains broad-band photometry of SN 1987A for 
t < 3 270 days. However, nearly all of the emission comes 
out in the far infrared which is not included in the ob- 
served bands. KF99 show that a more definite estimate 
can be made if one could measure the flux in a few iron 
lines, mainly [Fe II] 25.99/xm. There are several reasons 
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Table 1. Log of observations. 



Date 


Phase" 


Instrument /Mode 


Exposure time'' 


96/05/26 


3 380 


SWS/AOT2 


7307 


96/06/08 


3 393 


SWS/AOT2 


2 545 


97/07/19 


3 799 


SWS/AOT2 


7307 


98/02/04 


3 999 


LWS/AOTl 


3 641 


98/02/04 


3 999 


SWS/AOTl 


6 750 



" Epochs in days past explosion. 
'' Time in seconds. 



for this. First, at ~ 4 000 days KF99 find that almost 
half of the luminosity from the supernova is emerging in 
this line. Second, because the iron lines are only little af- 
fected by freeze-out they are good tracers of the instanta- 
neous energy deposition, and the flux of [Fe II] 25.99^m 
is therefore almost proportional to M(^'*Ti). Third, this 
line is formed through coUisional excitations and is there- 
fore insensitive to uncertainties in atomic data involved in 
calculating the recombination cascade. In a preliminary 
analysis, Borkowski ct al. (1997) find that their Infrared 
Space Observatory (ISO; Kessler et al. 1996) observations 
give an upper limit on the line flux which corresponds to 
only 1.5 X IQ-^ M©. 

Here we report on observations we have made with 
ISO. Although we have observed the supernova over the 
entire wavelength range 2.38 — 197/xm, wc will concen- 
trate our discussion on the region, 23 — 27/xm, where 
the strongest emission lines are expected to emerge. We 
compare these observations with theoretical modeling. Wc 
also briefly discuss our observations at longer wavelengths 
where we detect emission originating from elsewhere in 
the LMC. 

2. Observations and Results 

We have used ISO to observe SN 1987A on several occa- 
sions. Both the Short Wavelength Spectrograph (SWS; de 
Graauw ct al. 1996) and the Long Wavelength Spectro- 
graph (LWS; Clegg et al. 1996) were used, and Table 1 
summarizes our observations. We will concentrate here 
mainly on the SWS observation from 4 February, 1998 
(Sect. 2.1), though we provide a consistency check of this 
observation against our other SWS observations. The LWS 
data are discussed in Sect. 2.2. 

2.1. SWS Observations 

The SWS observation on 4 February, 1998 was made in 
the SWSOl mode with speed 4 and was centered on the 
position of SN 1987A, i.e., R.A. = 5h 35m 28.G5s; Decl. = 
-69° 16' 11'.'64; (J2000.0). SWSOl provides spectra which 
together cover the entire wavelength region between 2.38— 
45.2/im. However, as our models predict [Fe I] 24.05/im 
and [Fe II] 25.99/xm to be by far the strongest emission 



ISO/SWSOl spectrum of SN 19B7A (day 3999). Bin size = 0.052 fim. 




23 24 25 26 

Wavelength {/^m) 



Fig. 1. ISO SWS/AOTl Band 3 spectrum of SN 1987A 
on day 3999. The spectrum was reduced using interactive 
software (see text). The bin size was set to 0.052/im, in 
accordance with the instrumental resolution. Overlayed 
(dashed line) is the modeled line emission from model 
M2 (without photoionization) described in Table 2 and 
Sect. 3.2. To obtain the modeled line profiles we have 
assumed a maximum core velocity of 2 000 km s~^, and 
that the emission is homogeneously distributed through- 
out the core (see Sect. 3.3). The modeled peak flux of 
[Fe II] 25.99/xm corresponds to what is needed for a 3a 
detection. 

lines from the supernova at this epoch (day 3 999; see Sect. 
3), we have concentrated on measuring the flux in the 
wavelength range including these two lines. This range, 
between 22.5 — 27.5/xm, is covered by band 3D of SWS. 

The reductions were made using the SWS Interactive 
Analysis software system (SIA) available at ISO Spec- 
trometer Data Center (ISOSDC) at the Max Planck Insti- 
tut fiir extraterrestrische Physik in Garching (MPE). The 
most recent set of calibration files equivalent to off-line 
processing (OLP/pipeline) version 7.0 was used. The in- 
teractive reduction allows special care to be given to dark 
subtractions, which is of particular interest when measur- 
ing low fiux levels. Flat fielding was also applied, but we 
have not made any fringe corrections since the fringes at 
low flux levels disappear in the noise. 

In Fig. 1 we present a fully reduced spectrum of SN 
1987A for band 3D. Although the nominal instrumental 
resolution of SWS is i? « 1 000, the slowest SWSOl mode 
with speed 4 degrades this by a factor of 2. We have there- 
fore averaged the spectrum with a bin size of 0.052/um, cor- 
responding to 600—650 km s~^ for the two lines of interest. 
As the lines could extend to well above ±2 000 km s^^ this 
resolution should be suSicient to resolve the lines. How- 
ever, neither [Fe I] 24.05/im nor [Fe II] 25.99/im are seen in 
the spectrum. The 'features' that do appear around ~ 23.5 
and 25.7/Ltm are most likely due to instrumental effects, as 
they are rather robust in the sense that they appear in 
many detectors and in both up and down scans. They are 



p. Lundqvist et al.: ISO SWS/LWS observations of SN 1987A 



3 



certainly not due to fringes, nor do we believe they are 
effects of pure noise. The zero level of the spectrum is well 
defined and a simple zero-order fit to the spectrum gives 
an RMS of 0.47 ,Jy over the range 22.5 - 27.5/im. 

Looking more closely at the spectral regions of the 
two lines, we have simply measured the RMS between 
23 — 25/im and 25 — 27/im separately. Again, wc used a 
zero-order fit to the smoothed spectrum. A zero-order fit 
should provide the most conservative way to estimate the 
'noise' in our data. The resulting RMS values are ^ 0.34 Jy 
and ~ 0.46 Jy for [Fe I] 24.05/im and [Fe II] 25.99/im, re- 
spectively. This gives Scr upper limits of ~ 1.02 Jy and 
^ 1.38 Jy for the peak of the profile of the two lines. 
These values do not change much even if we double the 
bin size. 

SN 1987A was also observed in the SWS02 mode on 
three occasions (see Table 1). This mode provides scans 
over shorter wavelength ranges than the full SWSOl scans. 
We have reduced also the SWS02 data using SIA. While 
[Fe II] 25.99/xm was looked at on all three occasions (see 
Table 1), [Fe I] 24.05/im was only looked at on June 8, 

1996. The flux calibration for the SWS02 mode is in gen- 
eral more reliable than for SWSOl at a specific wavelength 
since dark exposures are taken immediately before and 
after the short SWS02 scans, and the sampling for the 
individual data points is better. However, the SWS02 ob- 
servations in Table 1 cover only the range 25.75— 26. 20/um, 
which could be just a fraction of the real line width. As the 
continuum level is therefore unknown, we have to rely on 
absolute flux calibrations, which are uncertain at the low 
flux levels of our data. Concentrating on the 26/im line, 
we measure ~ 0.75 Jy and ^ 0.62 Jy for the two epochs 
with the longest exposures, May 26, 1996 and July 19, 

1997, respectively. The accuracy of absolute fluxes for this 
wavelength regime at higher flux levels is < 20% (ISOSWS 
Data Users Manual, v5.0), but for our very low flux levels, 
~ 50% is probably more fair. 

Due to this uncertainty we can only say that our 
SWS02 flux estimates are consistent with the upper limits 
we obtain from the SWSOl measurement. 

2.2. LWS Observations 

Our LWS observations were made on 4 February, 1998. 
The full range of the LWS was scanned, i.e., 45 — 197/um. 
Standard reduction of the spectra using the ISO Spectral 
Analysis Package (ISAP) software resulted in the com- 
bined spectrum displayed in Fig. 2. Just as for the SWS 
observations, no emission was observed from the siipcr- 
nova. However, four unresolved emission lines from neigh- 
boring photoexcited gas were identifled: [O I] 63.2/im, 
[O HI] 51.8/im, [O III] 88.4/xm and [C II] 158Aim, with the 
fluxes (1.5±0.3) X IQ-^^ erg s"! cm'^, (1.39±0.13) x 10"" 



ISO/LWS spectrum of SN 1987A (day 3999). 
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lines sit on top of a broad continuum which peaks at 



100 120 140 

Wavelength (yum) 

Fig. 2. ISO LWS spectrum of SN 1987A on day 3999. 
The continuum up to 140/im can be fitted with a spec- 
trum emitted by dust with a temperature of ^ 37 K (solid 
smooth line). Longward of 140/zm a second component ap- 
pears to add in. The temperature characterizing this com- 
ponent is ^ 10 K. The combined emission is shown by the 
dashed line. See Sect. 2.2 and 4.3 for details. 



~ lOOjUm at a level of (17 ± 5) Jy, and which is only 
marginally weaker in the range 145 — 165//m. The in- 
tegrated continuum flux is (6 ± 2) x 10" erg s^^ cm^^, 
which for a distance of 50 kpc corresponds to a luminosity 
of (1.8 ±0.6) x 10^* erg. It should be pointed out that the 
spatial region sampled has a large diameter (~ 80"), which 
at a distance of 50 kpc distance corresponds to ~ 19 pc. 
We discuss the LWS results briefly in Sect. 4.3. 



3. Interpreting the SWS observations 

To understand what the upper limits on the iron lines in 
Sect. 2.1 mean, we have made model calculations. The 
models are described in Sect. 3.1, and we discuss the re- 
sults from the calculations in Sect. 3.2 



3.1. The model for line emission from the supernova 

The code we have used for our theoretical calculations is 
described in detail in Kozma & Fransson (1998a), and here 
we will only give a brief summary of the model. 

The thermal and ionization balances are solved time- 
dcpendcntly, as are also the level populations of the most 
important ions. 

The ions included are H I-II, He I-III, C I-III, N I-II, O 
I-III, Ne I-II, Na I-II, Mg I-III, Si I-III, S I-III, Ar I-II, Ca 
I-III, Fe I-V, Co II, Ni I-II. The following ions are treated 
as multilevel atoms: H I, He I, O I, Mg I, Si I, Ca II, Fe 
I, Fe II, Fe HI, Fe IV, Co II, Ni I, Ni II. The remaining 
lines are solved as 2-level atoms or as 4 — 6 level systems. 
For Fe I we include 121 levels, Fe II 191 levels, Fe III 112 
levels and for Fe IV 45 levels. A total of approximately 
6 400 lines are included in the calculations. 
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The radioactive isoptopes included are "''^Ni, "'''Ni, and 
44rpj^ These radioactive decays provide the energy source 
for the ejecta, and we calculate the energy deposition of 
gamma-rays and positrons solving the Spencer- Fano equa- 
tion (see Kozma & Fransson [1992] for a more detailed 
description of our thermalization of the gamma-rays). We 
have assumed that the positrons deposit their energy lo- 
cally, within the regions containing the newly synthesized 
iron. 

We have included 0.07 M© of ^^Ni (Sect. 1) and 3.33 x 
10~^ Mq of ^"^Ni, which corresponds to a ^''Fe/^^Fe ratio 
equal to two times the solar ratio. This mass of ^''Ni gives 
a good fit to the bolomctric light curve when the effects 
of freeze-out are taken into account (Fransson & Kozma 
1993). It is also consistent with observations of [Fe II] and 
[Co II] IR hncs (Varani ct al. 1990; Danzigcr ct al. 1991), 
observations of the ^''Co 122 keV line (Kurfess et al. 1992), 
and theoretical nucleosynthesis calculations (Woosley & 
Hoffman 1991). In our 'standard' model we have used l.Ox 
10-4 Mq of "^"Ti (e.g., Kumagai et al. 1991; Timmes et 
al. 1996). We have also made calculations with M{'^'^Tij = 
0.5 X 10-4 Mq and M{'^'^Ti) = 2.0 x 10-" Mq. 

Until recently, the lifetime of "^^Ti has been very poorly 
known. In our calculations we have used an e-folding time 
of 78 years thought to be a representative value to those 
found from experiments. However, while completing our 
calculations, we were informed about the recent measure- 
ments mentioned in Sect. 1 (S. Nagataki, private commu- 
nication). The new, and accurate value, 87.0 ±1.9 years, 
is close to what we have used, but decreases the instanta- 
neous radioactive energy input from the ^4^^^ (j^j^ay by 
^ 12%. This systematic error, shifting our estimated 
value of M(44Ti) upward by a similar amount, has been 
included in the error analysis in Sect. 4.1. 

As input for the calculations we have adopted the 
abundances from the lOH explosion model (Woosley & 
Weaver 1986; Woosley 1988). Spherically symmetric ge- 
ometry has been assumed, with shells containing the ma- 
jor composition regions as given by the explosion model. 
We have used the same distribution of shells as in Kozma 
& Fransson (1998a) where hydrogen is mixed into the 
core. The iron-rich core is assumed to extend out to 
2 000 km s^^, outside of which we attach a hydrogen en- 
velope, out to 6 000 km s~^. 

For the line transfer we have used the Sobolev approx- 
imation. This is justified for wcU-scparatcd lines in an ex- 
panding medium. However, especially in the UV, there are 
many overlapping lines, and one can expect UV-scattering 
to be important. The importance of line scattering de- 
creases with time as the optical depths decrease. The effect 
of line scattering is to alter the emergent UV-spectrum, 
but it also affects the UV-ficld within the ejecta. The ion- 
ization of elements with low ionization potential is sensi- 
tive the UV-field. During scattering the UV photons are 
shifted towards longer wavelengths, both due to a pure 
Doppler shift, but also because of the increased probability 



of splitting the UV-photons into several photons of longer 
wavelengths. A more accurate treatment of the line scat- 
tering is therefore expected to decrease the importance of 
photoionization. For the modeling of [Fc II] 25.99/im, the 
photoionization of Fe I to Fe II is likely to introduce the 
largest uncertainty. To check the importance of this effect 
we have in some models simply switched off the photoion- 
ization. A more comprehensive line transfer modeling is 
discussed in KF99. 

3.2. Model calculations 

The results of our calculations can be seen in Table 2. We 
tabulate line fluxes for four strong lines for three values 
of M(44Ti): 5 x IQ-^ (Ml), 1 x IQ-* (M2) and 2 x 10-^ 
Mq (M3). For each value of M('*''Ti) wc show results for 
models with and without photoionization. In all six models 
a simple form of dust absorption was assumed (Kozma 
& Fransson 1998b). The line fluxes in Table 2 are for a 
box line profile between ±2 000 km s-^, and assuming a 
distance to the supernova of 50 kpc. 

From Table 2 it is evident that [Fc II] 25.99/im is by 
far the strongest line in our simulations, and that its flux 
scales roughly linearly with MC^^Ti), as expected from 
Sect. 1. From power-law fits to the results of the calcula- 
tions described in Table 2, we obtain /26Mm oc M{'^'^Tiy-^^ 
when photoionization is included, and nearly the same 
without photoionization (/26/im oc M('''*Ti)^-^^). The sec- 
ond strongest line is [Fe I] 24.05/im. Its dependence on 
M(44Ti) is weaker: /24^m oc Mi^^Tif-^*^ and /24^^ oc 
^(44rjj^o.73^ respectively. Other lines from SN 1987A are 
far too weak for our ISO observations. 

The stronger dependence of /26/im on MC^^Ti) than 
for /24/im is due to ionization: while a higher M('*'*Ti) 
boosts the relative fraction of Fe II, Xpeii, the relative 
fraction of Fe I, Xp^i, decreases. For example, in models 
with photoionization, the average value of Xpcii in the Fe- 
rich gas increases from ~ 0.32 to ~ 0.55 when M{*'^Ti) is 
increased from 5 x IQ-^ to 2 x IQ-*, while Xp^i decreases 
from ^ 0.65 to ~ 0.36. The fact that the 26/im line is 
so much stronger than the 24/L(m line despite the roughly 
similar relative fractions of Fe I and Fe II is just a result 
of the much larger collision strength of the Fe II line than 
for [Fe I] 24.05/xm. 

Table 2 shows that switching off photoionization does 
not have a dramatic effect on the line fluxes. In the case 
of M2, /26/im decreases by ^ 15% and /24/im increases by 
28% when we turn off photoionization, simply due to a 
minor shift in ionization between Fe I and Fc II. As wc will 
discuss in Sect. 4.1, there are other uncertainties which are 
of the same magnitiude. 

3.3. Mass of^^Ti from [Fe II] 25.99iJ.m 

With our results it is straightforward to estimate the up- 
per limit on M(4'*Ti). Using f^&nm ^ 1-38 Jy from Sect. 
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Table 2. Modeled line flux in Jy at 4000 days"'*" 



Line 




M2'* 


MS'' 


[Fcl] 24.047//m 
[Fell] 25.995^im 
[Fell] 35.359Mm 
[Felll] 22.930Mm 


0.11 (0.12) 
0.45 (0.40) 
0.01 (0.01) 
0.01 (0.01) 


0.16 (0.20) 
1.03 (0.88) 
0.04 (0.04) 
0.03 (0.03) 


0.23 (0.33) 
2.19 (1.89) 
0.12 (0.11) 
0.11 (0.11) 



" Distance = 50 kpc. Box line profile of width ±2 000 km s 

' For values in parantheses photoionization has been excluded. 
" M(^*Ti) = 5 X 10-5 Mq 
M(**Ti) = 1 X 10"* Mo 
" M{*^Ti) = 2 X 10"" M© 



2.1, and our power-law fit to the results in Table 2 (Sect. 
3.2), we find an upper limit which is M{^Ti) < 1.3(1.5) x 
IQ-^ M0 with (without) photoionization included. These 
masses are for a box-shaped line profile, and therefore 
provide conservative limits. However, there is good rea- 
son to believe that the line profile should be similar to 
those for the lines observed by Haas et al. (1990) ~ 400 
days after the explosion. The strongest line observed by 
Haas et al. was [Fc H] 17.94/im and it had wfwhm = 
(2 900 ± 80) km s~^. Although it peaked just short off 
+1000 km s"^, and thus was not symmetric around the 
rest velocity of the supernova, it had the general ap- 
pearance of the expected line profile formed by a filled 
sphere with constant emission throughout the sphere. For 
a sphere extending out to 2 000 km s"-'^, the peak is 1.5 
times higher than for the box profile used in Table 2 (which 
is valid for a hollow sphere), and Ufwhm = 2 828 km s~^. 
Using such a line profile in our upper limits on M(*^Ti), 
we instead obtain M(44Ti) < 0.9(1.0) x lO"'' M© with 
(without) photoionization included. Within the frame- 
work of our modeling, a conservative limit (i.e., the case 
when photoionization is unimportant) for a plausible line 
profile is therefore M^^Ti) < 1.0 x lO""^ Mq. In Fig. 1, we 
have included the expected line emission for such a model 
with this limiting mass of '*'*Ti. We will evaluate this limit 
on M(44Ti) in Sect. 4.1. 

4. Discussion 

4-1- Uncertainties in the modeling 

There are several uncertainties involved in our modeling 

of the line fluxes. We have already checked the effect of the 
lifetime of ^^Ti (Sect. 3.1). We have also studied the effect 
of photoionization, and found that it introduces rather 
mild uncertainties. A similar level of uncertainty is due 
to the distance to the supernova. This is still inaccurate 
to the level of 5 — 10% (Lundqvist & Sonneborn 1999; 
Walker 1999), which means an uncertainty in the line flux 
of ~ 10 - 20%. 



Atomic data of iron arc notoriously difficult to calcu- 
late accurately. This is therefore another source of error 
in our modeling, especially for individual lines. Most im- 
portant for the 26/im line is the collision strength of that 
transition, r226/im- One normally assigns an uncertainty in 
the collision strength for the strongest iron lines to ~ 30% 
(e.g., Verner et al. 1999). In our models we have used 
f^26/im = 5.8 (Zhang & Pradhan 1995). To study the effect 
in detail we have tested a model with M(*^Ti) = 10"* M© 
and il26tim = 2.9. We find that there is not a linear scaling 
between /26/im and ri26/im, as one might naively believe. 
Instead, we find from linear interpolation that a 30% de- 
crease in collision strength gives a ^ 17% lower /26/jm, and 
a correspondingly higher estimate of M(^^Ti). The reason 
for this is that the gas is slightly hotter in the model with 
reduced collision strength, boosting the exponential term 
in the coUisional rate so that it somewhat counteracts the 
reduced collision strength. We have recently found out (A. 
Pradhan, private communication) that the preferred value 
for r226/im at the low temperatures in SN 1987A is most 
likely closer to ^ 7.0, i.e., higher than we have used. De- 
spite this, we have generously assigned an uncertainty of 
15% (upwards) in M(^^Ti) due to atomic data. 

Another source of uncertainty could be the explosion 
model used. In these calculations we have used the abun- 
dances from the lOH explosion model. A comparison be- 
tween the two models lOH (Woosley & Weaver 1986; 
Woosley 1988) and llEl (Shigeyama et al. 1988) was done 
in Kozma & Fransson (1998b). There it was found that the 
iron lines are not sensitive to the explosion model used, 
because the iron core mass is the same in both models. 
The iron core mass is set by the amount of ^®Ni which 
is accurately determined from the bolometric light curve. 
The choice of explosion model thus does not seem to be 
a major source of uncertainty when modeling these iron 
lines. 

In our calculations we assume a local deposition of 
the positrons originating from the radioactive decays of 
M(*^Ti). We believe this is a good approximation since 
optical and near-IR light curves of Fe I and Fe II lines show 
that trapping must occur (Chugai et al. 1997; Kozma & 
Fransson 1998b). Actually, there is no obvious sign of a 
leakage of positrons, neither from broad-band lightcurves 
(KF99), nor from the optical Fe I lines at 6 300 A until the 
last data point at 3 597 days in Kozma (1999). Although 
the trapping may well be fully complete, we have assigned 
an error to this assumption by 15%. 

Another approximation in our models is the assump- 
tion of a homogeneous density in each Fc-rich shell of the 
model core. To test the sensitivity to this assumption, we 
have run a model similar to M2, with photoionization in- 
cluded (sec Table 2), but where we have divided the mass 
in the Fe-rich ejecta into two components of equal mass 
but with different densities. The denser component is set 
to be nearly five times more dense than the other. De- 
spite the significantly different density distribution in this 
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model compared to that in M2, the differences in /26/im 
and /24/im between the models are small. Instead of the 
values listed in Table 2, the fluxes in the two-component 
model are 0.98 and 0.20 Jy, rcspccively, i.e., /26/^m differs 
by only ~ 5% compared to that in M2. This indicates that 
the model assumption of a homogeneous density in each 
Fe-rich shell docs not introduce a major uncertainty. 

Finally, screening and cooling by dust are potential 
sources of error in our models. The effects of dust are 
examined in Kozma & Fransson (1998ab). The screening 
we use is discussed in Kozma & Fransson (1998b), and is 
based on estimates by Lucy et al. (1991) and Wooden 
et al. (1993). Dust from pure iron is unlikely to form, 
as that would cool and block out all iron line emission 
once the dust has formed. On the contrary, there is a 
wealth of iron lines from the core at late times. KF99 
and Kozma (1999) find good agreement between modeled 
and observed broad-band lightcurves for t ^ 3 270 days 
and spectra at 2 870 days. As was pointed out also for 
the positron leakage, the models are also able to repro- 
duce optical Fe I lines at even later epochs. Based on this, 
we believe dust effects are small enough to neglect in our 
estimate of M('"^Ti). 

None of the model approximations we have used ap- 
pears to be uncertain enough to allow /24;^m to be of the 
same magnitude as f^^^m- The best estimate of M(^'*Ti) 
(within the framework of our modeling) should therefore 
come solely from f26fj.m- To estimate the combined error 
of f26nm due to all model approximations (except for the 
dust distribution in the ejecta), we adopt the uncertain- 
tics 15%, 15%, 15% and 5% for photoionization, distance, 
atomica data and clumping, respectively. For the choice 
of input model and positron leakage we adopt 15% each. 
This gives a combined uncertainty which is ^ 34%. On 
top of this we must add the maximum systematic error of 
12% discussed in Sect. 3.1 for the lifetime of *^Ti. With 
the line profile discussed in Sect. 3.3, we therefore arrive 
at an upper limit on M(''''Ti) which is ~ 1.5 x 10~^ M©. 
We note that this limit excludes the upper ends of the 
allowed ranges of M{^'^Ti) found by Chugai et al. (1997) 
and KF99 (see Sect. 1). Combining our limit with the pre- 
liminary results of KF99 for the broad-band photometry, 
a likely range for M{^^T\) is (0.5 - 1.5) x lO"'' M©. We 
emphasis, however, that the lower limit of this range is 
probably more uncertain than the upper (for the reasons 
mentioned in Sect. 1), which is indeed indicated by the 
preliminary analysis of Borkowski et al. (1997). 

4-2. Implications of the derived mass of^Ti 

Models for the yield of *^Ti give quite different results. 

This is most likely due to how the explosion is generated 
in the models, and how fallback onto the neutron star is 
treated. Timmes et al. (1996; see also Woosley & Weaver 
1995) use a piston to generate the explosion, and they ac- 
count for fallback in a rather self-consistent way. In their 



model with zero-age mass Mzams = 20 Mq (i.e., corre- 
sponding to SN 1987A) the mass of the initially ejected 
^^Ti is 1.2 X 10-4 M0, but of this only 1.4 x 10"^ Mq es- 
capes after fallback. This is less than we argue for in Sect. 
4.1, and could suggest that fallback was not important for 
SN 1987A, though we caution again that the lower limit 
found by KF99 (see also Kozma 1999) may not be very 
strict. If fallback is unimportant the ejected amounts of 
^^Ni and ^^Ni would be too high in this model, typically 
by a factors of ^ 2 — 4, judging from the effects of fall- 
back in the 25 model in Woosley & Weaver (1995). It 
should be emphasized that the variation of M{'^'T\) with 
.^ZAMS in Timmes et al. (1996) is complex, and that for 
models with Mzams = 18 M© and Mzams = 22 Mq, the 
calculated M{^T\) comes within the range we propose, 
albeit close to our lower limit. 

The models of Thielemann et al. (1996) simulate the 
explosion by depositing thermal energy in the core, and 
they insert the mass cut artificially so that the right 
amount of ejected ^^Ni is produced. (This effectively 
means that fallback is included also in these models.) Sim- 
ulating the explosion in this way ensures larger entropy 
and thus more alpha-rich freeze-out than in Woosley & 
Weaver (1995). Accordingly, the ratio M {'^'^Ti) / M {^^m) 
(where M (^^Ni) is the mass of ejected ^^Ni that docs not 
fall back) is higher in the models of Thielemann et al. than 
in piston-driven simulations. For example, in the 20 M© 
model of Thielemann et al. (1996) M(^^m) « 0.074 Mq, 
M(5^Ni) « 2.9 X 10^3 a,nd M{^^Ti) « 1.7 x lO""^ Mq, 
with M(^'^Ni) defined in the same way as M(^^Ni) and 
M(44Ti). The values of M{^^m) and M{^'^m) are close 
to what have been inferred for SN 1987A (Suntzeff & 
Bouchet 1990; Fransson & Kozma 1993). The titanium 
mass is slightly larger than the upper limit of the range we 
estimate in Sect. 4.1. So, while our estimate of M{^^T\) 
cannot rule out with certainty any of the two methods 
used for the explosion (piston-driven or heat generated), 
our results could indicate that an intermediate method 
should be used (see also the discussion on this in Timmes 
et al. 1996). From models of the chemical evolution of 
the Galaxy, and especially the solar abundance of ^^Ca, 
a value for M{^^T\) closer to that of Thielemann et al. 
(1996) may be more correct, at least for supernovae in 
general. 

In this context we note that a higher value of Mi^'^Ti) 

is produced in asymmetric explosions (Nagataki et al. 
1998). The method of calculation employed by Nagataki 
et al. (see Nagataki et al. 1997) is similar to that in Thiele- 
mann et al. (1996), though the models of Nagataki et al. 
allow for 2-D instead of just 1-D. With no asymmetry, the 
models of Nagataki et al. produce M{^T\) ~ 6 x 10"^ Mq 
for an explosion similar to SN 1987A, when the mass cut 
has been trimmed to M(^^Ni) « 0.07 Mq. This is signifi- 
cantly less than Thielemann et al. (1996) despite the simi- 
lar method of modeling. Applying an asymmetry by a fac- 
tor of 2 between the equator and the poles, the explosion 
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energy becomes concentrated toward the equator resulting 
in relatively stronger alpha-rich freeze-out, which increases 
M(44Ti) to - 2 X 10--* Mq (for M{^^m) « 0.07 Mq). 
This is outside our observationally determined range and 
could indicate that asymmetry was not extreme in SN 
1987A, especially if M{'^'^Ti) is close to the limit derived 
by Borkowski et al. (1997). We note that a piston-driven 
calculation could perhaps allow for asymmetry, as such 
models give very small values of M(4'*Ti) in 1-D. 

A direct way to estimate M(*'^Ti) in supernovae is to 
observe the gamma-ray emission from the radioactive de- 
cay of ^^Ti. The 1.156 MeV hue associated with the decay 
of '^^Ti has only been detected in two supernova remnants 
(and no supernovae): Cas A (lyudin et al. 1994; The et al. 
1996) and the newly discovered gamma-ray/X-ray source 
GRO J0852-4262/RX J0852. 0-4622 (lyudin et al. 1998; 
Aschenbach 1998). While M (^**Ti) in the latter is not yet 
well-known, MC^^Ti) in Cas A was 1.7t^j x 10"^ Mq 
(Gorres et al. 1998). In the models of Nagataki et al. 
(1998), the value of M('*''Ti) at the upper end of this 
range would indicate that Cas A exploded asymetrically. 
Our estimate of MC'^Ti) in Sect. 4.1 does not seem to in- 
dicate a high degree of asymmetry for SN 1987A, but an 
independent measurement of its 1.156 MeV line may be 
needed to be more conclusive on this point. To detect this 
emission from SN 1987A, instruments hke INTEGRAL 
(Leising 1994) are required. 

4-3. Interpretation of LWS observations 

The two [O III] lines observed (see Sect. 2.2 and Fig. 2) 
can be used to estimate the average density of the emitting 
gas. We have used a multi-level atom to do this, using the 
atomic data of Mendoza (1983), Osterbrock (1989) and 
Aggarwal (1993). Assuming a temperature of 10"* K, we 
obtain a mean electron density of 120 ± 75 cm^'^ in the 
[O III] emitting gas. 

The continuum can be fitted with a spectrum of the 
form Fx oc Bx{T){l — cxp{—Tu)) , where td is the dust op- 
tical depth. With a functional form of td — ro.55(0.55/A)" 
(where A is in microns) we obtain a best fit with T ~ 37 K, 
ro.55 ~ 0.74 and a ~ 0.83. The fit (see Fig. 2) works well 
up to ~ 140/im, where an extra source appears to add in. 
This could indicate the presence of cold dust. Assuming 
the same functional form for this emission as for the 37 
K component, the temperature of the cold component is 
close to 10 K. 

5. Conclusions 

We have looked at SN 1987A with ISO SWS/LWS. The 
supernova was not detected in any of the spectra. In 
particular, we have derived upper limits for the fluxes 
of [Fe I] 24.05/Ltm and [Fe II] 25.99^m, and made time 
dependent calculations for the late line emission from 
the supernova. We have assessed various uncertainties to 



the models, and for a plausible line profile we then find 
an upper limit on the mass of ejected '^''Ti, M{*'^Ti) ~ 
1.5 X 10^"* Mq. Together with the preliminary results of 
KF99 this brackets the ejected mass of '^^Ti to the range 
M(44Ti) = (0.5 - 1.5) X 10-4 Mq, which is close to the 
yield of '''^Ti in models of the explosion by Woosley & 
Weaver (1995), Thielemann et al. (1996) and Nagataki et 
al. (1997). The lower limit is probably less stringent than 
the upper, as indicated by the results of Borkowski et al. 
(1997). A more direct hmit on M(4''Ti) can be made when 
instruments measuring the gamma-ray line emission from 
the "'''Ti decay will be available. 

The only emission our ISO observations detect is from 
gas and dust in the direction toward the supernova. In 
particular, the [O III] 51.8//m and [O III] 88.4/im lines 
indicate a gas density of 120 ± 75 cm"'^, and the dust 
continuum can be explained by a temperature of ~ 37 
K. A second dust component with ^ 10 K may also be 
present. 
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